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ABSTRACT 

"a^" 

Qs I We present results from detailed imaging of the centrally dominant radio elliptical 

galaxy in the cooling flow cluster Abell 2597, using data obtained with the Wide 
Ch ' Field and Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). This 

object is one of the archetypal "blue-lobed" cooling flow radio elliptical galaxies, also 
displaying a luminous emission-line nebula, a compact radio source, and a significant 
dust lane and evidence of molecular gas in its center. We show that the radio source 
is surrounded by a complex network of emission-line filaments, some of which display 
^ i a close spatial association with the outer boundary of the radio lobes. We present 

fT^ \ a detailed analysis of the physical properties of ionized and neutral gas associated 

J^ I with the radio lobes, and show that their properties are strongly suggestive of direct 

(^ ' interactions between the radio plasma and ambient gas. We resolve the blue continuum 

O^ ■ emission into a series of knots and clumps, and present evidence that these are most 

■^ ■ likely due to regions of recent star formation. We investigate several possible triggering 

Oj. mechanisms for the star formation, including direct interactions with the radio source, 

X. , filaments condensing from the cooling flow, or the result of an interaction with a 

|H \ gas-rich galaxy, which may also have been responsible for fueling the active nucleus. We 



(N 



propose that the properties of the source are plausibly explained in terms of accretion 
of gas by the cD during an interaction with a gas-rich galaxy, which combined with the 
fact that this object is located at the center of a dense, high-pressure ICM can account 



XI 

5_j ■ for the high rates of star formation and the strong conflnement of the radio source. 
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1. Introduction 



Rich clusters of galaxies generally contain a diffuse intracluster medium (ICM) having 
temperatures T ^ 10^ — 10^ K and densities n ~ 10~^ 



10 ^ cm ^, as inferred from observations 



at X-ray and other wavebands (e.g., Sarazin 1986; Fabian 1994). An important question related 
to the overall physics of the ICM concerns the central regions of clusters (r ^ 10 — 100 kpc), 
where the inferred ICM densities and pressures in some cases are sufficiently high that cooling to 



T <^ 10 K can occur on timescales shorter than the cluster lifetime (e.g., Cowie & Binney 1977; 



Fabian & Nulsen 1977). These "cooling flow" clusters often exhibit intense optical emission- line 
nebulae associated with the centrally dominant (cD) galaxies at their centers, together with blue 
continuum excess emission, and the strength of these effects appears to correlate with the cooling 
rate or central pressure of the X-ray emitting gas ( Peckman 1981 ; Hu, Cowie, & Wang 1985; 
Johnstone, Fabian, Sz Nulsen 1987 ; Romanishin 1987| Romanishin fc Hintzen 1988 ; Heckman 
et al. 19891 ; [McNamara fc O'Conneh 1989| , |1992| , |1993| ; [Crawford fc Fabian 199J , |1993| ; [Allen 1995 ; 



McNamara 1997 ). The majority of cooling flow cluster cD galaxies are also radio-loud, hosting 
moderately powerful but relatively compact radio sources ( Burns 1990| ). 



Several intriguing questions have been posed concerning cDs in cooling flows. For example, 
one question concerns the energy source for the line emission: the optical line luminosities are 
generally one to three orders of magnitude higher than would be expected from the X-ray-derived 
cooling rates, if each proton in the cooling gas experienced only a single recombination (e.g.. 



Heckman et al. 198£; Voit & Donahue 1997). This has led to the consideration of numerous other 



potential sources of heating, including photoionization by young stars or the AGN, auto-ionizing 



shocks, photoionization by X-rays from the ICM, or cosmic-ray heating (e.g., Voit & Donahue 



1997[ and references therein). Likewise the blue continuum emission is of considerable interest, and 



suggested possibilities for its origin have included scattered nuclear light, recent star formation, or 
inverse Compton emission. Finally, a crucial point to investigate concerns the inter-relationships 
between these various components: the cD elliptical, the emission-line gas, the blue continuum, 
the radio source, and the surrounding dense high-pressure ICM: how are these related to one 
another, and do they have underlying common physical relationships? 

We present results from a detailed study of a comparatively nearby strong cooling flow cluster, 
Abell 2597, using images obtained with the Wide Field and Planetary Camera 2 (WFPC2) on 
board the Hubble Space Telescope (HST). Previous radio and X-ray images have revealed that the 
cluster contains a moderately luminous radio source (-Pi. 4 gHz ~ 6 x 10^^ W Hz~^, pwen, Ledlow 
fc Keel 1995[; Zem = 0.0821 ± 0.0002, [Voit fc Donahue 1997t ^abs = 0.0823 ± 0.0003, [Taylor et al." 



1999 ) embedded in a luminous X-ray halo with Lx ~ 10^^ ergs s ^ ( prawford et al. 1989| ; pwen 



White, fc Burns 1992 ; Sarazin et al. 1995| ; Sarazin fc McNamara 1997[ ). The host galaxy of the 



radio source is a cD elliptical containing an extended, luminous emission-line nebula (Heckman 



let al. 19891 ; [Smith fc Heckman 1989| ; [Crawford fc Fabian 1992t [Voit &: Donahue 1997[ ). The 



galaxy is remarkable for its "blue lobes" of continuum excess emission that are approximately 
co-incident with the radio lobes ( [McNamara &: O'Connell 1992[ , |1993| ). Spatially extended 21-cm 
H I absorption has also been detected against the radio lobes ( P'Dea et al. 199^ . Although the 
radio source is comparatively compact (with a size ^ 5 kpc), its power is comparable to that of 
much larger FR II ( "classical double" ) sources, suggesting that the dense intracluster medium may 
play a role in its confinement. The source may be the nearest example of a compact steep-spectrum 
source ( p'Dea et al. 1994 ). Here we present data that fully resolve the blue emission in the lobes 



and reveal a complex network of emission-line and continuum filaments surrounding the central 
source. 

All distance-dependent quantities in this paper have been calculated assuming 
Hq = 75 km s~^ Mpc"^ and qq = 0. 



2. Observations 

We observed A2597 in July 1996 with the HST WFPC2 Camera #3 (WF3) through the 
F410M filter, with the total exposure time being "CR-split" into three equal, consecutive exposures 
to facilitate the removal of cosmic rays. This is effectively a narrow-band image, mostly containing 
flux from the redshifted [O II] A3727 emission line. In this paper we also discuss new work that we 
have carried out on two archival broad-band images of A2597, taken with the WFPC2 Planetary 
Camera (PCI) through the F450W and F702W filters (Program ID #6228 in the HST Data 



Archive; see also Holtzman et al. 1996 ). These images contain substantial continuum flux in 
addition to the redshifted [O II] A3727 and Ha + [N II] AA6548,6583 emission lines, and therefore 
allow a detailed investigation of the emission-line and stellar continuum properties of this object, 
when combined with the F410M image. In Table || we present further details of the observations 
and filter characteristics. 

Recalibration of the data was performed using the IRAF/STSDAS[| pipeline processing 
software for WFPC2 (see also Holtzman et al. 1995a, [1995b[ ), using updated calibration reference 



files where applicable. The removal of hot pixels with elevated dark currents was investigated 
using the IRAF/STSDAS task WARMPIX together with a table of hot pixels nearest in time to 
the observations, although it was found that most of these were adequately removed through the 
use of the appropriate dark frame. 

We verified that no significant alignment shifts were present between the multiple exposures 
in each of the three filters, before combining them using the IRAF/STSDAS task CRREJ which 
removes most of the cosmic rays as well as some remaining hot and cold pixels. The cosmic ray 



'^The Image Reduction and Analysis Facility (TRAP) is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc., under contract 
to the National Science Foundation. The Space Telescope Science Data Analysis System (STSDAS) is distributed by 
the Space Telescope Science Institute. 



and bad pixel removal was checked by examining the residual difference between the cleaned 
image and each of the two original frames, and interpolations were performed across the remaining 
defects. Sky subtraction in the broad-band F450W PCI image was carried out by obtaining a 
mean value of the sky background in several empty regions of the frame, and subtracting this 
as a constant value from the entire image. The sky background for the F702W PCI image was 
determined by scaling the F450W sky value using the sky count rates presented for these filters in 
the HST Data Handbook. This was necessary since the extended envelope of the elliptical is more 
prominent in the F702W image than in the F450W image, thereby precluding an accurate direct 
determination of the sky background. The sky count rate values (0.0071 and 0.0028 e~ s~^ pixel"^ 
for the F450W and F702W images respectively) were sufficiently low that residual large-scale 
fluctuations in the images were at levels ^ 0.05 counts, well within the Poisson noise limits. 

The two broad-band images were transformed to F450W and F702W magnitudes, m,F45ow 
and mpjQ2W respectively, in the HST WFPC2 synthetic "VEGAMAG" photometric system, 
defined in Poltzman et al. (1995b|) , making use of the IRAF/STSDAS SYNPHOT package to 
obtain accurate photometric parameters for the filters used. Conversions to Johnson B and R 
magnitudes were not carried out at this point, since substantial differences between the WFPC2 
filter response functions and those of the Johnson standard system would introduce photometric 
errors if significant color variations are present across the galaxy. 



3. Morphological Properties 

In Figure |I] we present the broad-band F702W PCI image of A2597. It is clear that the 
central region of the elliptical is dominated by a complex network of filaments. Some of the 
filaments appear clumpy and marginally resolved, particularly the arcs extending toward the 
north-east and the south. There is also some evidence of a dust lane crossing the nucleus, at 



Table 1. HST WFPC2 Imaging Observations of A2597 



Date 



Exposure 
(s) 



Filter 



(A) 



5A^ 

(A) 



Detector Pixel size 
(arcsec) 



1996 July 27 


2100 


F410M 


4090.1 


146.7 


WF3 


0'.'0996 


1995 July 5 


2500 


F450W 


4519.1 


956.6 


PCI 


0'.'0455 


1995 July 5 


2100 


F702W 


6867.8 


1382.5 


PCI 


0'.'0455 



^The specifications for the mean central filter wavelength A and filter width 
6X are given in Table 6.1 of the WFPC2 Instrument Handbook ([Biretta et al 



1996 ), where definitions for these quantities are also presented. 
Pixel sizes are from Holtzman et al. (1995a| ). 



approximately the same position angle as the major axis of the elliptical. A small object, possibly 
a companion galaxy, is present to the north of the elliptical, and numerous star clusters are also 
distributed throughout the elliptical, some of which may be resolved (cf. [Holtzman et al. 1996| ). 

The filamentary structure dominates the narrow-band emission-line image of the galaxy, as 
shown in our F410M WF3 image presented in Figure H which contains the redshifted [O II] A3727 
emission. We superpose contours from a 3.5 cm VLA image of the radio continuum emission 



( Sarazin et al. 1995 ). The radio and optical images were aligned on the basis of the radio core 
position, which is assumed to coincide with the broad-band optical nucleus of the galaxy. Having 
aligned the images on this basis alone, it is remarkable that the bright northern filamentary 
arc follows precisely along the outer contours of the radio lobe. The rest of the northern radio 
continuum is also apparently associated with lower-level filamentary emission. Although there 
is less optical filamentation associated with the southern radio lobe, an emission-line filament is 
present near the eastern edge of this lobe. Closer to the center, there is a bright, compact region of 
emission immediately south-west of the nucleus, which is co-incident with a radio hotspot as well 
as an abrupt change in the orientation of the radio lobe. We also note that some of the filamentary 
structures within the galaxy are apparently unassociated with the radio lobes, particularly the 
filaments extending away from the nucleus toward the east, and a more diffuse region surrounding 
the nucleus on the west. 



3.1. Separating the Continuum and Emission-Line Features 

One of the interesting questions related to the properties of the filaments concerns the relative 
fractions of continuum and emission-line contributions to their total flux. We have investigated 
this by comparing the narrow-band F410M image to the broad-band F450W image, both of which 
contain the redshifted [O II] A3727 emission line but contain different amounts of stellar continuum 
flux. We used the IRAF/STSDAS package DRIZZLE ( |Fruchter fc Hook 1997| ) to transform the 



higher-resolution PCI F450W image to the coarser pixel scale of the WF3 F410M image, also 
taking into account the differences in orientation angles and geometrical distortions across the 
chips. 

We used the SYNPHOT package to flux-calibrate each of these two images in cgs units 
(ergs s~^ cm~^ arcsec^^). A fundamental uncertainty inherent in this calibration is the lack of a 
detailed knowledge of the spectral energy distribution across the filter bandpass at each pixel in 
the image. Two limiting cases are: (1) monochromatic emission, with all the flux being emitted 
at one wavelength (appropriate for pixels dominated by a single emission line); (2) a flat spectral 
energy distribution, speciflcally with constant flux density per unit wavelength fx. For each of 
the two fllters, we generated a flux calibration factor for each limiting case using the SYNPHOT 
task BANDPAR. The difference between the two cases is negligible for the F410M filter (~ 1.4%) 
but is of order ~ 12% for the broader F450W filter. Therefore we flux-calibrated the F410M 
filter using the monochromatic limit, while the F450W filter was calibrated with the flat-spectrum 
approximation, which is justified by the fact that the flux in this broader filter is dominated by 



stellar continuum (based on examining the spatially integrated ground-based spectra of Voit & 



Donahue 1997). We note however that this approximation may break down for individual pixels 
with extremely high [O II] A3727 equivalent widths, introducing a maximum error of ~ 12% (for 
pure line emission), or more likely ^ 5 — 6% (for comparable proportions of line and continuum 
emission). 

We then subtracted the flux-calibrated F410M image from the geometrically transformed 
F450W image to investigate whether the filaments display any residual continuum emission. As 
described above, the uncertainties in fiux calibration due to differences in the filter bandpasses 
are likely to be ^ 5%, therefore this subtraction provides a reliable indication of the relative 
contribution of emission-line and continuum emission to the filaments. We found that the extended 
filaments and arcs around the radio lobes disappeared almost entirely in the subtracted image, and 
therefore can be considered to consist predominantly of [O II] A3727 line emission. The filaments 
nearer the galaxy nucleus, however, display a mixture of continuum and line emission and are 
apparently dominated in some regions by excess blue continuum, which we discuss in more detail 
in the following section. 



3.2. The Blue Continuum Excess 

We have used the broad-band PCI F450W and F702W images to investigate the detailed 
spatial distribution and morphology of excess blue emission in A2597, which was first discovered in 
ground-based images by McNamara &: O'Connell (1993| ) as marginally resolved U — I continuum 



excess aligned along the radio axis. We first matched both PCI images to a common point-spread 
function (PSF) by convolving each image with the PSF of the other, having computed the 
appropriate PSF for each filter band using the HST Tiny Tim software ( Krist fc Hook 1997| ). 



The images were then calibrated in the HST VEGAMAG system, after which we calculated the 
"T'F450W — ''7^F702W color distribution of the galaxy, which we present in Figure ^ We have also 
verified that smoothing the HST data to the resolution of the ground-based images reveals a good 
correspondence between the f7-band and F450W blue continuum. The blue emission in Figure ^ 
is resolved in detail into a number of separate components: (1) a band of extended, intense blue 
emission around the "base" of the northern radio lobe, immediately north-east of the nucleus; 
(2) a group of blue point-like features south-west from the southern radio lobe; (3) a blue arc that 
appears to trace the edge of the southern radio lobe; (4) diffuse blue emission around the radio 
lobes, and specifically to the north-west of the nucleus. 

The possibility that the blue continuum is predominantly due to scattered nuclear light is 



ruled out on the basis of its low observed polarization ( [McNamara et al. 1998D ; optical synchrotron 
or inverse Compton emission are also not feasible since none of the blue continuum sources 
correspond to strong radio hotspots. The most plausible remaining explanation is that the blue 
continuum is due to young stars, and in this section we examine various star formation scenarios 
that may be responsible for producing the extended blue continuum and the isolated blue clusters. 
In §P we present a more detailed discussion of the relative plausibility of several possible triggering 
mechanisms for the star formation, together with various scenarios for the origin of the gas from 
which the stars are formed. 



3.2.1. The Extended and Diffuse Blue Continuum 

Our HST images resolve an intense band of blue continuum to the east and north of the 
nucleus, extending around toward the west, thereby covering a spatial scale of several kiloparsecs 
(Figure |^. Some less intense blue emission is also associated with the southern radio lobe. 
Resolving the blue emission allows a more accurate estimate of its intrinsic color and luminosity 



than was previously obtainable from the ground-based measurements obtained by McNamara & 



O'Connell (1993 ), therefore we use our WFPC2 data to investigate the star formation processes 



involved. 

The intrinsic color of the blue excess is derived from the residual continuum in each band 



after subtracting the underlying elliptical fit, as described in more detail in §3.3. In that section 
we also discuss the effects of dust obscuration, which is present in the F702W band as well as 
F450W. Since the strongest obscuration is confined to a few small localized patches, we are still 
able to discuss the overall properties of the remainder of the blue continuum excess. Both the 
F450W and F702W bands are also contaminated by redshifted emission lines, [O II] A3727 and 
Ha+[N II] respectively. However, the total integrated flux from the excess blue emission in the 
F450W image is substantially higher than the flux from the [O II] A3727 emission in that region 
observed through the F410M filter; this also applies to the comparison between the F702W excess 
continuum and the Ha+[N II] flux scaled from the [O II] A3727 flux using the mean Ha+[N II] / 
[O II] A3727 ratio for this object. Furthermore, the bulk of the line emission displays little spatial 
correspondence with the blue continuum. Therefore in this analysis we neglect to first order the 
contribution from line emission to the continuum color; spatially resolved spectrophotometry 
would be required to further quantify its contribution. We measure the intrinsic color of the blue 
excess to be in the range r?T,F45ow — "iF702W ~ 0-4 — 0.8 with a surface brightness ranging from 
"1F450W ~ 22.0 — 20.8 mag arcsec"^; the integrated absolute magnitude of the entire extended 
region of blue continuum excess is Mf45ow ~ —18.7 (spread over a projected area of ~ 20 kpc^). 

We compare the measured properties of the blue continuum to a set of stellar population 



synthesis models ( Chariot &: Bruzual 199S , private communication), which have been K-corrected 
to the redshift of A2597 and cahbrated in the HST VEGAMAG system using our HST WFPC2 
filter bandpass transmission curves. We find it preferable to retain our measured magnitudes 
in the VEGAMAG system and compare them with models that have also been calculated in 
this system, rather than attempting to transform the magnitudes to other systems such as the 
Johnson/Cousins system, because uncertainties in the color terms can translate into errors of 
several tenths of a magnitude. The models are shown in Figure ^ and represent two scenarios 
for star formation: (1) a constant star formation rate (SFR), assuming an invariant IMF; (2) an 
instantaneous burst of star formation (or "Simple Stellar Population", SSP). In each case we 



investigated two standard IMFs for which models were available, namely [Salpeter (1955 ) and 



Scalo (19861) , each with upper and lower mass cutoffs of 0.1 and 100 Mq respectively. We also 



examined models of different metallicities, specifically Zq, 0.2Zq and 0.02Zq. Further details 



of the model parameterizations are described in Chariot & Bruzual (1991, 1999) and Bruzual 



fc Chariot (1993 ). Having calculated HST colors and magnitudes for these various models, we 
note that the difference between the color evolution of Zq and 0.2Zq models is minimal for the 
constant SFR models, thus in Figure ^ we only plot the 0.2Zq models for the purposes of clarity. 



For the SSP models, the color evolution is more sensitive to metallicity but is relatively insensitive 
to the choice of IMF, therefore the SSP models that we plot are for the Salpeter IMF only. 

In Table ^ we present the inferred star formation parameters that are required to match 
the observed range of colors and integrated absolute magnitude of the extended blue continuum 
distribution, for the various star formation scenarios that we have described. For each model, we 
tabulate the upper and lower limits of the range of epochs in its evolution that reproduce our 
observed mF45ow — "if702W color. For constant SFR models, the SFR for each of the two epochs 
is obtained by matching the model to the Mf45ow absolute magnitude. For the instantaneous 
burst models we tabulate the total mass of stars (an SFR is not physically meaningful in this case 
since the burst is represented by a delta- function) . 

We point out that the inferred ages are all highly dependent upon the measured colors, and 
should therefore be interpreted as upper limits due to the likely presence of reddening across most 
of the field. For example, a decrease in the intrinsic mF450W — ''7^F702W color by only a further 
0.2 magnitudes will yield ages as low as 6 x 10^ yr for the continuous SFR models, and ~ 5 x 10^ yr 
for the instantaneous burst models. However, due to our lack of accurate reddening estimates 
(precluded by the presence of emission lines in both bands) we are unable to quantify this further. 

The southern radio lobe displays substantially less blue emission than the northern lobe. Its 
W'F450W ~ 'n'F702W color is similar to that of the northern lobe but its flux is much less, with an 
integrated apparent magnitude of only ?7iF450W ~ 23.2. Since it is fainter, errors in its measured 
properties are much more dominated by uncertainties in the underlying continuum, reddening 
correction, and contributions from line emission, therefore we do not feel justified in carrying out a 
comparison with star formation models to the same level of detail as the more prominent northern 
blue excess. We do note, however, one interesting feature of the southern blue emission, namely a 
blue "arc" that follows the outer contours of the southern lobe, from its western edge around to 
the south. This is unlikely to correspond to optical synchrotron emission, as the radio contours 
show no enhanced brightening at the location of the arc, thus the arc may also indicate emission 
from young stars although our data are of insufficient S/N to quantify this further. 



3.2.2. The Isolated Blue Knots 

We now discuss a group of four blue clumps, or knots, that are readily apparent south-west of 
the southern radio lobe. Because of their relative isolation, they can serve as a valuable reference 
for comparison when discussing the star formation that has occurred north of the nucleus. In the 
original (unsmoothed) PCI images they are marginally resolved, with FWHM values in the range 
~ Gf!l2 — C/.'IG. They apparently correspond to the southern U — I blue continuum excess reported 
by [McNamara fc O'Connell (1993| ). These features do not appear to be directly associated with 
the radio lobe, and in all likelihood they are young star clusters. In Table ^ we present the 
magnitudes and colors of these objects, calibrated in the HST VEGAMAG photometric system. 
Our photometric analysis of these objects indicates that their colors are about 1.5 — 2 magnitudes 
bluer than the surrounding elliptical continuum. As a consistency check, we also carried out 
photometry of a number of the brighter, redder objects around A2597 in the F702W image, and 



compared these with the i?-band photometry of the same objects presented by Holtzman et al. 
(1996). We found that the magnitudes typically agreed to within ~ 0.1 — 0.2 mag, which is quite 



satisfactory given the uncertainties in the color terms that had been used by Holtzman et al. 
(1996) to transform their HST magnitudes to the Johnson/Cousins system. 



We have compared the colors and magnitudes of the blue clumps to the stellar synthesis 
models of pharlot fc Bruzual (1999| ), K-corrected to the redshift of A2597 and calibrated in 
the HST bandpasses as described in the previous section. Since we are investigating individual 
clusters in this case, we carry out a comparison only with the instantaneous burst models, with 
the same range of metallicities and IMFs as for the northern blue continuum. In the lower panel of 
Figure ^ we plot the evolution in color-magnitude space of three population synthesis models with 
a Salpeter IMF, metallicities of Zq, Q.2Zq and 0.02Zq, and total mass 10^ -^0; together with 
the observed colors and magnitudes of the blue clusters. The comparison suggests that the colors 



Table 2. Star Formation Parameters for the Extended Blue Continuum 



Model 




log(T) (yr) 


SFR (M0/yr) 


Const. SFR: Salpeter: 


Zq 


8.8-9.7 


1.1-0.8 




0.2Zq 


8.7-9.7 


1.0-0.6 




O.O2Z0 


9.2-10.1 


0.7-0.4 


Scalo: 


% 


8.6-9.3 


2.1-1.0 




0.2^0 


8.6 - 9.4 


1.7-0.7 




O.O2Z0 


8.9-9.8 


1.0-0.4 


log(M/M0) 


Inst. Burst: Salpeter: 


Zq 


7.6 - 8.7 


8.5-9.3 




0.2Zq 


7.4-8.8 


8.3-9.3 




O.O2Z0 


8.3-9.0 


8.8-9.1 


Scalo: 


Zq 


7.7-8.7 


8.7-9.3 




0.2Zq 


7.5-8.9 


8.6-9.3 




O.O2Z0 


8.4-9.1 


8.9-9.3 



Note. — The upper and lower limits of the timescale r of 
each model correspond to the epochs in the model evolution 
that produce the observed range of mF45ow — ''tT'F702W colors. 
For the Constant SFR models, the tabulated SFR is that which 
is required to match the observed integrated magnitude of the 
excess blue continuum. For the Instantaneous Burst models, we 
instead tabulate the total mass of stars since in this model all the 
star formation occurs at one instant in time. Further details are 
given in the text. 
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of the three bluest clusters are consistent with ages ^5 — 10 Myr and nietallicities above ~ 0.2Zq, 
with ages of the two redder of these three clusters ranging up to a maximum of ~ 100 Myr for 
masses up to ~ 5 x 10^ -^0- The reddest of the four clusters is either much older (~ 10^ yr) 
with a mass ~ 10^ -^0! or otherwise is the same age as the other three, but has higher extinction 
and/or a metallicity significantly above solar (however, this cluster is fainter than the other three 
and its properties are not as well constrained). In § ^ we discuss in further detail the implications 
of the measured properties of the blue continuum regions for various physical mechanisms that 
can be responsible for triggering the star formation. 



3.3. Elliptical Isophote Fits 

In order to present a more detailed discussion of the various components that make up the 
central galaxy, including the filaments, the blue continuum, and possible dust obscuration, it 
is necessary to model and remove the underlying elliptical stellar continuum distribution. This 
was done by means of the IRAF/STSDAS task ELLIPSE (based on an algorithm described by 



Jedrzejewski 1987), which makes use of iterative sigma-clipping algorithms to reject the possible 



contribution of non-elliptical components to the fit. 

Our data are complicated both by the presence of filaments in emission, and also by apparent 
dust absorption in the central regions of the galaxy. We initially attempted techniques similar 
to those described by Carollo et al. (1997) and [Forbes, Franx, fc Illingworth (1995| ), involving a 



number of repeated fits to both the F702W and F450W images, which allows the possibility of 
correcting for patches of dust obscuration. However, this technique does not easily allow additional 



Table 3. Compact Blue Objects in A2597 

ID Aa A6 mF450W "T-F702W WF450W - "IF702W 

(") (") (mag) (mag) (mag) 



1 


-1.79 


-2.21 


26.41 ± 0.48 


25.55 ±0.26 


0.86 ±0.55 


2 


-3.39 


-1.94 


25.38 ±0.16 


25.07 ±0.14 


0.31 ±0.21 


3 


-2.14 


-2.99 


25.42 ±0.18 


25.10 ±0.16 


0.32 ±0.24 


4 


-3.47 


-2.57 


25.64 ±0.23 


25.88 ±0.29 


-0.24 ±0.37 



Note. — For each object, the offsets in Right Ascension and 
Declination (Aa and A6 respectively) are given in arcseconds relative 
to the continuum nucleus of the galaxy. The 1 a uncertainties on the 
individual magnitudes r?T,F45ow and 7tt,f702W are derived from the count- 
rate statistics, and are combined in quadrature to yield uncertainties for 
the color index mF45ow — "^F702W • 
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corrections for emission features, particularly the line emission from the filaments which is present 
in both the F450W and the F702W broad-band images of A2597. 

Therefore, we obtained the underlying elliptical continuum by fitting isophotes to a 
near-infrared H-band image of A2597, recently obtained with the F160W filter on NICMOS by 
Donahue et al. (1998). This band has the advantage of containing no line emission from the 



filaments, and also with the effect of dust obscuration being much reduced (although still present). 
We fitted elliptical isophotes to this image by allowing the flux, ellipticity and position angle of 
each isophote to vary freely, and using sigma-clipping to reject features deviating by more than 3 a 
from the mean level in each isophote (for example corresponding to stars or clusters, or the dust 
obscuration region). Inspecting the resulting residual image showed that this model provided a 
good fit to the overall galaxy profile, with the only significant residuals corresponding to localized 
features such as clusters or dust absorption patches. Some low-level residual structure on large 
scales was also present, which is attributed to uncertainties in the fiat-fielding applied in the 
reduction process. 

The elliptical model obtained from the NICMOS data was then applied to the F702W and 
F450W images, by keeping the isophote shape parameters constant (i.e., ellipticity and position 
angle), but allowing the fiux of each isophote to vary freely, and using sigma-clipping as before to 
remove contributions due to the filaments or dust patches. For isophotes with a semi-major axis 
^ 6" we had to allow the isophote shape parameters to vary, since the smaller size of the NICMOS 
image, together with its low-level fiat-field effects, produced inaccuracies in the model at these 
large radii. However, for isophotes inside this radius the NICMOS isophote shapes provided a very 
good fit to the F702W and F450W data. Subtracting the resulting model from the input image 
yielded for each filter a first-iteration residual image which contained most (but not yet all) of the 
"excess" emission. Because of the substantial contributions in emission from the filaments within 
the central ~ 1 — 2", we carried out a second iteration, according to the following procedure. 
Pixels deviating from by more than 3 times the rms value in the residual image were subtracted 
from the original input image to yield a second input image. New isophotes were then fitted to 
this image, using the isophote shape parameters from the first iteration and allowing only the fiux 
to vary, and again using sigma-clipping to reject non-elliptical emission features. This yielded a 
second residual image; repeating this procedure further showed that the process had converged by 
the second iteration to provide an accurate representation of the underlying elliptical continuum 
flux distribution, together with the residual emission and absorption features. 

In Figure |^ we present the final set of isophotal parameters fitted to the galaxy. The elliptical 
continuum color index across most of the galaxy is in the range m,F45ow — "^F702W ~ 1.5 — 2 
which corresponds to 5— F~l — 1.5,a typical value for gE galaxies (e.g., parollo et al. 1997| ). 
The color is somewhat redder in the central region, possibly related to the presence of dust. One 
interesting feature of the galaxy is a change of ~ 90° in the orientation of the major axis, relative 
to the outer isophotes, for semi-major axis lengths ^ 1" . This also corresponds to an increase in 
ellipticity toward the center. This change may be due to the presence of dust located on either side 
of the nucleus, although in the F160W image the change is quite smooth and therefore we cannot 
rule out the possibility that this may refiect a genuine shift in the underlying light distribution, 
for example due to a nuclear bar (e.g., Nieto et al. 1992|) . 
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In Figure ^ we present the "residual" filamentary structure in the F702W band, produced 
by subtracting the elliptical galaxy model from the original image. We also overlay contours 
representing the VLA 3.5 cm radio continuum image ( Sarazin et al. 1995| ). The positive 
emission features in this image are likely to be dominated by Ha line emission rather than stellar 
continuum, based on the fact that the filamentary structures in the F410M image (Figure |2|) are 
dominated by [O II] A3727 emission. The filaments in Figure g are sampled with both higher 
S/N and better spatial resolution than those in Figure |2|: the bright arcs are resolved into clumps 
in several locations, and diffuse emission can be clearly distinguished across the face of the 
northern radio lobe. In addition, the presence of strong, patchy dust obscuration is evident in a 
number of locations, particularly to the south-east and to the north of the nucleus. Some weaker 
obscuration is also present west of the nucleus, just beyond the bright radio hotspot and associated 
emission-line filaments. 



3.4. The Emission-Line Morphology 

The distribution of optical line emission around the northern lobe of the radio source, as 
displayed in Figures ^ and ^ can be divided into three regions: (1) the bright "upper" arc along 
the north-west of the lobe; (2) a less distinct "lower" arc around the southern boundary of the 
lobe; (3) diffuse, lower surface-brightness emission across the face of the radio lobe. The upper 
arc runs parallel to a bright ridge of enhanced radio emission along the edge of the lobe, both of 
which display a linear extent of ~ 2" (2.8 kpc). The arcs are dominated by emission-line radiation; 
the lack of any detected blue continuum emission excludes non-thermal processes such as optical 
synchrotron emission or inverse Compton scattering. We discuss these features in more detail in 



Table 4. Emission-Line Fluxes and Luminosities 



Region 



Fr 



[O II] A3727 



(ergs s 



'^ cm ^1 



L[0 II] A3727 

(ergs s~^) 



Lr/3 

(ergs s~ 



N. arc - upper 
N. arc - lower 
N. lobe face 
N. lobe (total) 
SW knot 
SE arc 



9.31 X 10"^^ 
9.49 X 10"^^ 
1.58 X 10"^'' 
3.46 X 10"^'' 
9.52 X 10-15 
8.45 X 10-15 



1.30 X 10^1 
1.33 X 10^1 
2.21 X 10^1 
4.84 X 10^1 
1.33 X lO^^i 
1.18 X 10^1 



2.96 X 10^° 
3.02 X 10^° 

5.02 X 10^° 
1.10 X 10^1 

3.03 X 10^° 
2.68 X 10^° 



Note. — The [O II] A3727 line fiuxes were obtained by 
summing all the pixels in the corresponding regions on the WF3 
F410M image, which is dominated by [O II] A3 72 7 line emission. 
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The southern lobe displays a compact knot of high surface-brightness line emission co-incident 
with the strong radio hotspot immediately south-west from the nucleus. A prominent semi-circular 
filament of emission is also present SE from the nucleus, along the eastern edge of the southern 
radio lobe. We find that these features are also dominated by line emission and not continuum 
radiation. We discuss the southern emission-line features further in §H. 

In Table |^ we present the integrated [O II] A3727 fluxes of each of the three northern 
emission-line regions, together with the total luminosity of emission-line gas surrounding the 
northern radio lobe, as well as the bright emission-line knot immediately SW from the nucleus 
and the semi-circular SE arc. We also derive the H/3 luminosities for each of these regions from 
the measured [O II] A3727 line fluxes by adopting a canonical value of H/3 / [O II] A3727 = 0.23 
(uncorrected for reddening), determined from the ground-based spectra of Voit &: Donahue (1997] ). 



Implicit in this conversion is the assumption that the excitation and extinction properties of the 
nebula as a whole are also representative of the material associated with each radio lobe; this 
approximation is necessitated by the current absence of spectra with higher spatial resolution. 



4. Interaction between the Northern Radio Lobe and Ambient Gas 

A prominent feature of this galaxy is the complex network of emission-line filaments associated 
with the radio lobes. In this section we examine scenarios in which the northern filaments trace 
physical interactions between the northern radio lobe and ambient gas. We derive physical 
parameters of the ionized and neutral gas around the radio lobe, including the neutral fraction 
and volume filling factor of the gas and its total mass, by combining constraints obtained from our 
new emission-line data with those from previous observations. 

We consider a geometry in which the radio lobe is surrounded by an approximately uniform 
shell of line-emitting gas (except along the edge furthest away from the nucleus, where no 
significant emission is detected). This scenario is suggested by the similarities in flux and spatial 
extent of the upper and lower edge-brightened arcs. Furthermore, the integrated flux of the diffuse 
emission across the face of the radio lobe (which would include emission both from behind and in 
front of the lobe) is approximately twice that of the flux in each arc, supporting an edge-brightened 
geometry for the shell. 

Limits on the shell thickness or depth d, obtained from the width of the arcs, are thus 
applicable to the line-emitting gas covering the face of the radio lobe. The width of the bright 
upper arc is only marginally resolved, covering less than ~ 2 PCI pixels (130 pc), therefore we 
retain d as a free parameter normalized to a fiducial value of 100 pc. The surface area of diffuse 
line emission across the face of the radio lobe is ~ 3.24 arcsec^, thus the volume of the emission-line 
region in front of the radio lobe is Vtot ~ 0.65 (d/lOO pc) kpc'^. Given standard assumptions such 
as Case B ionization conditions, uniform density emission-line gas, and an edge-brightened shell 
geometry, our observed H/3 luminosity then provides a combined constraint on the volume filling 
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factor /jj+ and electron density rie of the warm ionized gas (with electron temperature Te): 

/ rp ^0.87/ J x-l 

nlf^+r^n -L^] ^^—] cm-6 (1) 

eJii+ Vio^k; vioopcy ' ^ ^ 

We aim to derive the physical properties of the gas by combining this constraint with those 
obtained from H I observations. [O'Dea et al. (1994 ) observed extended, broad H I absorption 



(FWHM ~ 412 km s~^) against the northern radio lobe. The kinematic properties are similar 
to those of the optical line-emitting gas, thus the northern radio lobe would be surrounded by a 
collection of clouds containing both neutral and ionized gas. In this section we derive physical 
properties of the gas for the following two plausible scenarios: 

1. The H I is co-spatial with the ionized gas, with each H I cloud possessing an ionized 



'skin" responsible for the line emission (e.g., O'Dea et al. 1994 ). In this case the ionized and 



neutral gas are swept up by the expanding radio lobe into a shell, thereby compressing the gas and 
determining its kinematics, but not necessarily ionizing it (alternative ionization mechanisms are 
discussed by [Voit Sz Donahue 1997 ). For simplicity, in this scenario we consider the line-emitting 



gas around the northern lobe to be entirely photoionized by radiation originating from an external 
source (e.g., the active nucleus, hot stars, or radiation escaping from shocks in another region). 

2. The ionized and neutral gas are located in different regions, for example if the northern 
radio lobe expands into a previously existing ensemble of neutral clouds (tidal debris or material 
condensed out of the cooling flow), fully ionizing each cloud by means of fast shocks driven into 
the gas. In this case, the radio lobe would be surrounded by an inner region of mostly ionized gas 



and an outer region of neutral clouds (see for example the schematic presented in Bicknell, Dopita 



fc O'Dea 1997| ). The ionized gas would contain a mixture of cooling "post-shock" gas together 



with photoionized "precursor" gas that absorbs ionizing radiation emitted by the shocked gas. 

4.1. Scenarios for the Distribution of Neutral and Ionized Gas 

4.1.1. Co- Spatial Neutral and Ionized Gas 

In this scenario, the neutral and warm ionized components occupy the same region but are 
represented by different volume filling factors. Parameterizing the depth of the H I absorbing 
material by d thus converts the observed H I column density (A^h ~ 4.5 x 10^*^ (Tg/lOO K) cm^^, 
O'Dea et al. 199^ ) to a combined constraint on the H I density and volume filling factor: 



^1 



{ Ts \ { d ^ _3 



Pressure balance is considered to hold between the neutral gas (spin temperature Tg ~ 100 K), 
warm ionized gas (Tg ~ 10*^ K) and the hot ambient medium (T -^ lO'^ K). We define ^ = Mjjo/Mjj+ 
to represent the relative masses of neutral and warm ionized gas, and combine this with equations 
(|l]) and (0) to obtain the gas densities: 

rn \ -I / rp s 0.87 

J- s \ I J-e \ „,„— 3 



"»^ ~ ^-'^l Took J loM?' 
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rp \ —2 / 'T' \ 1-87 

"HO - 1.5 X 103 d :r7^ ) -r^] cm-3. (3) 



100 Ky Vio'^K 

Taking nH+ ~ ng -^ 200 cm"^ ( |Voit fc Donahue 1997| ) to represent the overall line-emitting gas, 
we obtain the volume filling factor and thus the total mass of neutral and ionized gas in the shell, 
assuming only that its thickness d remains approximately uniform: 

M„. = 3.7X10= (^)°" (^5^)" Mo. 

M„. = LOxlO«(^) Mq, (4) 

If the gas is distributed in the form of approximately spherical neutral clouds with ionized 
skins, then the above constraints yield the relative extent of the H"*" and H'' cloud radius: 






1 + 7.3 X ' '^ ^ ' "" ^ ' 



ViOOK; VlO^Ky V200cm 



3 

(5) 



We use this, together with the projected surface area of the northern lobe and the H I volume 
filling factor, to constrain the H I cloud covering fraction Cf and cloud radius: 



rno Cf = 2.7 xW-'(-^) (-^) f ,,,";i,3 ) pc . (6) 



-1 
VlOOKy VlO^Ky V200cm~ 
Combining this with the limits on the covering fraction of the extended H I absorption in front 



of the northern radio lobe {cj > 0.006, O'Dea et al. 1994) and using our determination of the 



surface area of line-emitting gas in front of the radio lobe, yields an upper limit on the cloud 
radius: rjjo < 0.46 pc , and corresponding numbers of clouds: N^i > 6.0 x 10^ . A covering fraction 
of unity is obtained with N^i ^ 2.8 x 10^^ clouds, with radii rjjo ^ 2.7 x 10^^ p^. 

We conclude our description of the physical properties of the gas in this 
scenario with a calculation of the inferred pressure associated with the gas: 
P ~ 5.5 x 10"^'^ (ne/200 cm"^) (Tg/lO'^ K) dyn cm~^ . This is in very good agreement 
with the central X-ray pressure Px ~ 6 x 10"^'^ dyn cm~^ derived by Sarazin et al. (1995| ) but 



is lower than the estimated minimum pressure of the northern lobe, Pme ~ 1.4 x 10~^ dyn cm~^ 
(re-calculated with the cosmological parameters used in the present paper), thereby suggesting 
that the radio lobe is somewhat overpressured with respect to the surrounding gas. 



4-1-2. Spatially Distinct Neutral and Ionized Components 

An alternative scenario for the distribution of neutral and ionized gas involves expansion 
of the radio lobe into a more extended region of neutral clouds. Each cloud is fully ionized by 
the passage of the bowshock, analagous to the scenario envisioned for GPS/CSS radio sources 
(Bicknell et al. 1997|). The shell of warm ionized gas around the radio lobe would be surrounded 



by neutral clouds whose origin may be unassociated with the radio lobe (for example, tidal 
debris around the host galaxy, or gas condensed from the cluster cooling flow). Line emission is 
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produced by cooling "post-shock" gas downstream from the shocks as well as "precursor" material 
immediately upstream from the shocks, photoionized by soft X-ray / UV radiation from the 
post-shock gas. The ionized gas clouds possess a bi-modal density distribution: the post-shock 
clouds have densities ^ 2 orders of magnitude above those of the photoionized precursor gas, 
while the relative total line radiation output of the two regions are approximately equal (Dopita 



fc Sutherland 19951) . 



We consider a case in which the expansion velocity of the emission-line shell surrounding the 
radio lobe is equated with a shock velocity V^h ~ 300 km s^^ (i.e., half of the observed optical 
linewidth). A relatively uniform expansion speed is suggested by the well-ordered morphology 
of the emission-line gas. MAPPINGS II models show that the post-shock material in this case 



cools to gas with a density ~ 230 times greater than the precursor material (Dopita & Sutherland 



1996; Koekemoer 1996; Koekemoer Sz Bicknell 1998). The electron density measured from the S II 



ratio by Voit &: Donahue (1997 ) represents an average over all the emission from both regions, 
thus we use their value together with the ratio of post-shock / precursor line emission to infer a 
precursor density ~ 2 cm~'^. Using our observed II/3 luminosity, together with the relationships for 
shock and precursor emission, yields a required shock surface area of ~ 25 kpc^ to account for the 
observed line emission. This compares well with the apparent surface area of the shell (~ 30 kpc^) 
— i.e., the covering fraction of shocks across the emission- line shell would be of order unity for 
precursor gas of this density. The corresponding spatial scale on which the precursor gas is ionized 
would be of order ^ 10 — 100 pc (depending on its volume filling factor); beyond this distance the 
gas would remain largely neutral. 



4.2. A Comparison of the Energy Budget for Each Scenario 

Having calculated the expected physical properties of the gas in each of the two preceding 
scenarios, we now compare their relative likelihood. In particular, we examine whether the radio 
lobe can supply the required energy to the surrounding gas. 

We begin with the total minimum energy content of the northern radio lobe, 
E^e = 1.95 X 10^^ ergs, the lower limit on the radio source age trad ^ 5 x 10^ yr from 
the synchrotron loss timescale, and the radio luminosity of the northern lobe '~ 4 x 10^^ ergs s~^ 
(derived from the observed radio properties by |Sarazin et al. 1995 , and rescaled to the cosmological 



parameters used in this paper). We assume a fiducial value of ^ 10% for the conversion efficiency 
of the total jet kinetic power Lj into radio emission, noting that this value is not well constrained 



observationally (|Eilek 1982t [Begelman, Blandford, fc Rees 1984| ), i.e., Lj ^4 x 10^^ ergs s ^. 



In the first scenario, the energy transfer from the expanding radio cocoon to the gas is 
predominantly kinetic. Given the cospatial nature of the gas components in this scenario, we 
equate the expansion speed of the optical line-emitting gas (300 km s~^) with that of the H I. The 
kinetic energy of the neutral and warm ionized gas components is then: 
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which is substantially less than the total minimum energy content of the northern radio lobe. 
In other words, the rate of energy input required to supply this amount of kinetic energy to the 
combined neutral and warm ionized gas components need only be ^ 5 — 10% of the total energy 
input into the radio lobe, therefore it is feasible that the expansion of the radio lobe can account 
for the gas kinematics of the northern emission-line filament system. 

For the second scenario, we investigate whether the energy input from the lobe expansion can 
feasibly account for the observed line luminosity through ionization by radiative shocks. If the 
rate of work done on the ambient medium by the expansion of the lobe corresponds directly to 
the total luminosity of the fully radiative shocks that are produced, then the resulting H/3 line 
luminosity can be written in terms of the jet energy flux and shock velocity (Bicknell et al. 1997): 



Lh« ~ 0.8 X 10^1 f P^^-^) °''Vi ^ A ergss-\ (8) 

"-"^ V300kms-V V4 X 10^3 ergs s-iy & ' ^' 

where we have equated the shock velocity Vsh with the 300 km s~^ expansion velocity of the 
emission-line shell surrounding the radio lobe. Despite the inherent uncertainties in deriving the 
radio lobe and jet energetic parameters, this value of Lh/3 is in good agreement with our measured 
total integrated line luminosity for the northern lobe, presented in Table 0. We note that this does 
not necessarily imply that shocks are the dominant ionization mechanism in the northern lobe; 
instead we are stating that shocks, if present, correspond to a set of physically plausible values of 
parameter space. There is some question whether the overall properties of the entire EELR around 
A2597 can be explained in terms of shocks, given the constraints on the electron temperature and 
gas ionization state presented by Voit &: Donahue (1997] ). However, their results apply primarily 



to the integrated spectrum of the emission-line nebula, to which the northern lobe contributes only 
~ 20% of the flux, thus local ionizing mechanisms in this region can readily be accommodated; 
higher spatial resolution spectra will prove useful in further addressing such questions. 

Thus, we have shown: (1) it is energetically feasible that the northern radio lobe can supply 
the kinetic energy of the surrounding gas, irrespective of the ionization mechanism; (2) if the 
radio lobe expansion indeed produces auto-ionizing shocks in the gas, then the required shock 
luminosity and other properties of the gas are in good agreement with the observations. 



5. Deflection of the Southern Radio Jet 

The bright radio hotspot O'.'S south-west of the nucleus is a region of great interest since this 
location corresponds to a sharp change in the direction of the radio jet toward the east, which also 
appears to dramatically affect the resultant morphology of the entire southern radio lobe, ^arazin 



et al. (1995 ) postulated the presence of a dense deflecting clump of cool gas (T ^ 10^ K), but the 
available (ground-based) optical data were of insufficient spatial resolution to verify its existence. 
Our HST images clearly resolve a knot of high surface-brightness optical line emission associated 
with the radio hotspot. This object is therefore a highly suitable candidate for conducting a study 
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of the effects of a direct jet / gas cloud interaction. Here we use our images to infer the physical 
properties of the gas clouds, including their size, mass and density (degree of inhomogeneity) . 

Two plausible scenarios exist for describing the encounter between the jet and the gas: either 
the radio jet propagates into a pre-existing region of cool dense gas, or an initially straight jet is 
deflected as a result of an interaction with one or more gas clouds moving directly into it (for 
example, due to their orbital motion around the galaxy). The number of free parameters in each 
case prohibits a discussion of which scenario is more likely. Rather, we examine fundamental 
physical properties of the gas that would be applicable to either scenario, using considerations 
based on energetics and momentum transfer, and following similar arguments to those previously 



presented by Eilek et al. (1984 ) and Sutherland, Bicknell, fc: Dopita (1993 ) 



We consider a non-relativistic jet with fiducial speed Vj ~ 10^ km s^^ and kinetic power 
Lj ~ 4 X lO'*^ ergs s~^. The acceleration imparted by the jet to gas approximately at rest is 
related to the jet momentum flux, thus the minimum gas mass required to signiflcantly deflect the 
jet over its lifetime (equated to the radio source age trad ~ 5 x 10^ yr) is M^ ^ 6 x 10^ -^0- If the 
gas cloud volume is approximately spherical (with a diameter given by its transverse dimension, 
del ~ 500 pc), the required mean density of hydrogen nuclei must exceed (n) ^ 4 cm~^. 

We next examine the heating of the gas by the jet, using en to represent the fraction 
of Lj which is converted into heat energy. The cloud needs to be dense enough to cool at 
this rate or faster, otherwise it would be heated to high temperatures and disperse. Writing 
Tirms as the rms average hydrogen number density of the cloud, and the cooling coefficient as 
A22 = A(r) / 10^^^ ergs s~^ cm^, the balance between the cooling and heating rates requires 

"™^"U22j UxlO«ergss-lj 1,500 pc J ' '•"' 

Comparing this with (n) suggests that the cloud either exceeds the lower limit on its mass and/or 
is inhomogeneous. Finally, the rms electron density ne,rms in the cloud that is required to produce 
the observed H/3 line luminosity is 

ne rms ~ H -^ —^ Cm'^ . (10) 

' v5oopcy vio^k; ^ ^ 

The similarity of the densities in equations (|9|) and ( p!o[ ) might suggest that the deflecting cloud 
was composed purely of warm ionized gas at 10'^ K. However, in photoionized gas at these 
temperatures the cooling coefficient is only A(T) ~ 3 x 10^^^ ergs s~^ cm^ and such gas could 
not radiate energy deposited by the jet unless the efficiency of deposition was low, en ~ 0.03. 
For typical conditions in photoionized emission line gas, the bolometric luminosity Lboi is about 
25 times that in H/3, or Lboi ~ 8 x 10"^^ ergs s~^, which is much smaller than Lj. Thus, either 
the efficiency of heating the deffecting cloud is low, or the cloud contains substantial amounts of 
atomic and/or molecular gas. The latter scenario is supported independently by the recent HST 
NICMOS detection ( Donahue et al. 1998| ) of signiflcant molecular H2 emission co-incident with 



the knot. Hence, we suggest that the observed properties of this emission-line knot are consistent 
with the presence of clumps of dense gas sufficiently massive to deflect the jet signiflcantly over 
the lifetime of the radio source. 
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6. Implications for Cooling Flows, cD Galaxies and Radio Sources 

We have so far demonstrated that the emission-line nebula in the cD elliptical of A2597 
displays a close association with the radio source, and that several different properties of the gas 
are strongly suggestive of direct interactions between the radio lobes and the ambient medium. 
We have also shown that the blue continuum emission originates in a number of resolved locations 
within the central few kpc, and does not display a direct spatial correspondence with the 
emission-line filaments. We have identified a number of regions of recent or on-going star formation 
that can be modeled either as single-burst populations with ages of several times 10^ — 10^ yr, 
or as continuous (constant) star formation with ages about an order of magnitude higher. We 
have also resolved an elongated region of strong dust obscuration immediately south-east of the 
nucleus, together with evidence of patches of dust in other locations nearby. 

In this section we explore the likely physical inter-relationships between these various 
components in cooling flow clusters, using A2597 as an example. Studies to date suggest that the 
strength of the blue continuum excess and nebular line emission around cD galaxies in cooling flow 



clusters increases as a function of the inferred total cooling rate ( McNamara fc O'Connell 1989 



Allen 1995 ; McNamara 1997] and references therein). Cooling flow cluster cDs possess abnormally 



high amounts of molecular gas and dust, and display more compact radio morphologies, relative 
to cD galaxies not associated with cooling flow clusters. Thus, it is crucial to understand the 
underlying physical connection between the dense, high-pressure ICM surrounding the elliptical, 
the blue continuum excess, the nebular line strength, the strong dust obscuration, and the radio 
source. In particular, is the presence of a cooling flow either a necessary and / or sufficient 
condition to explain these properties of the central galaxy, or do other processes also play a role 
(e.g., mergers or interactions)? We address this issue by breaking it down into the following 
questions: (1) what has triggered the recent star formation; (2) what is responsible for triggering 
the AGN; (3) how are the triggering mechanisms related to the fuel supply in each case, and what 
is the relationship between the star formation mechanism, the triggering of the radio source, and 
the presence of a cooling flow? 



6.1. What has Triggered the Star Formation? 

Generally, star formation in dense molecular clouds is triggered by compression of the clouds 
resulting from enhanced external pressure, or by deep molecular shocks associated with cloud-cloud 
collisions. Here we examine the likelihood of various means of inducing these effects to result in 
the formation of stars. 

One triggering mechanism for star formation is direct interactions between the radio jets or 



lobes and dense gas clouds immediately surrounding them (Rees 1989; Begelman & Cioffi 198S 



De Young 1989, 1995; Daly 1990| ). This scenario invokes the collapse of marginally unstable clouds 



in the ambient medium, resulting from overpressurization as they are overtaken by the expanding 
radio cocoon. It predicts that, while star formation is initiated along the radio lobe boundaries, 
the spatial correlation is short-lived due to differences between the orbits of the stars and the 
motion of the radio plasma. The relatively mild expansion velocity of the cocoon, compared to 
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the working surface of the jet, prevents disruption of the clouds and faciUtates their cohapse. The 
large-scale asymmetry of the blue continuum in A2597 would correspond to an intrinsic asymmetry 
in the initial distribution of molecular gas clouds. The radio source age in this scenario would have 
to agree with the ages of the young stars, which we find to be in the range ~ 10^ — 10^ yr from the 
instantaneous burst stellar model synthesis ages. The lower limit on the radio source lifetime is 
^ 5 X 10^ yr from the synchrotron spectral ageing analysis, but given the uncertainties involved in 
this derivation the actual radio age can be substantially higher. Thus, from timescale arguments 
it is plausible that the star formation is triggered by the radio source. This would imply that the 
radio source is intrinsically old but unable to penetrate the dense surrounding ICM to scales of 
more than a few kpc. 

Star formation in cooling flows is also proposed to occur via pressure-driven collapse of 



filaments or clouds condensing out of the cooling flow (Fabian et al. 1986; McNamara & O'Connell 



1989; Allen 1995). The inferred cooling flow mass accretion rates for A2597 (Crawford et al 



1989; Edge et al. 1992) are about two orders of magnitude above our values of the SFR. Such 



discrepancies are well-known from previous studies of A2597 and other cooling-flow clusters. 
Possible solutions include erroneous estimates of the SFR, for example due to a time- variable SFR 
or changing IMFs; or mass deposition into material that does not necessarily form stars. These 
and other possibilities are discussed in detail by Johnstone, Fabian, &: Nulsen (1987| ); |0'Connel] 
fc McNamara (19891) ; [Thomas fc Fabian (1990|) ; and [McNamara fc Jaffe (1994|) . We do not rule 
this out as a potential source of fuel for star formation. It does not readily predict, however, the 
highly asymmetric structure of the blue continuum in A2597. One means of obtaining asymmetry 
is via subcluster merging activity (e.g., prawford &: Fabian 1992 , [1993| paines, Fabian, fc Thomas 



1994). However, A2597 displays little evidence of a recent cluster/subcluster merger, as the 



X-ray data suggest a relaxed system ( parazin &: McNamara 1997[ ). The spatial scale of the blue 
continuum emission is also one to two orders of magnitude smaller, r ^ 3 kpc, compared with the 
typical cooling radius of ~ 50 — 150 kpc at which the stars are expected to form (Johnstone et al. 



1987 ). Thus, cooling of the X-ray gas can provide a source of cooler material for star formation, 



but an additional mechanism appears to be required for triggering the concentrated massive star 
formation observed in A2597. 

This leads us to discuss star formation triggered via accretion of a gas-rich companion galaxy 
by the elliptical. Galaxy interactions are generally associated with triggering star formation (e.g., 
[Toomre &: Toomre 1972 ; see also [Balick &: Heckman 1982 , and Barnes & Hernquist 1992 for 
reviews). Detailed n-body/hydrodynamic models show that interactions between ellipticals and 
small gas-rich galaxies generally lead to gas accretion into the elliptical core ( [Weil fc Hernquisi 
[1993 ), triggering a centrally condensed starburst as well as star formation in the tidal debris 
(e.g., Mihos & Hernquist 1994a, 1994b| , [1996 ). Such features are observed in a wide range of 
galaxies, but specifically in large ellipticals (e.g., ptrom 



Strom 1979; McNamara & O'Connell 



1992 ; Holtzman et al. 199(: ; Sparks, Carollo, &: Macchetto 1997| ). The star- forming regions in 



A2597 show similar characteristics: intense star formation near the galaxy core, with smaller, 
more isolated blue clusters further out. The primary objection to transporting gas-rich galaxies 
into the cluster core is that ram-pressure stripping by the hot ICM removes much of the relatively 
diffuse H I gas ( punn fc Gott 1972 ; Nulsen 1982| ). However, dense clouds of molecular gas and 
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dust have much lower filhng factors and higher column densities than the H I component and are 
less affected by hydrodynamical processes ( [Kritsuk 1983 ; Valluri fc Jog 1990| ). Indeed, studies of 



the molecular content of galaxies in clusters show no significant dependence on location in the 
cluster, whereas the H I content decreases substantially toward the cluster core (Kenney & Young 



1989; Casoli et al. 199l| [Lavezzi fc Dickey 1998 ). Therefore such an interaction can quite plausibly 



provide a fresh supply of gas, predominantly molecular. The stellar model ages of ~ 10^ — 10^ yr 
are also typical of the dynamical timescales involved in such interactions. 

Thus, we propose that two plausible triggering mechanisms for the star formation are: 
(1) related directly to the passage of the radio source, or (2) related to accretion of a gas-rich 
companion galaxy. If the star formation is modeled in terms of a single burst then the resulting 
ages are in the range ~ 10^ — 10^ yr — the low end of this range is comparable to the lower 
limit on the radio source age, thus we associate the single-burst models with triggering by the 
radio source. If, on the other hand, the star formation is modeled using a constant SFR, then the 
required ages of the stars are much higher, of order a few times 10^ — 10^ yr. This is too high to be 
accounted for in terms of radio source triggering, but is consistent with the dynamical timescales 
typically required for interactions between galaxies. While our current data do not allow us to 
distinguish further between these two possible triggering mechanisms, we are able to conclude that 
triggering by the radio source is consistent with a single burst ~ 10^ — 10^ yr ago, while triggering 
by an interaction corresponds to continuous, on-going star formation over longer timescales up to 
~ 10^ yr. 



6.2. AGN Trigger Mechanism, Fuel Supply, and Origin of the Gas Filaments 

Having discussed the star formation trigger mechanism, we now investigate the trigger for the 
AGN itself. The majority of cD ellipticals in cooling flow clusters are known to host radio-loud 
sources, compared to a much lower fraction of cD's in non cooling flow clusters ( [Burns 1990 ). A 



number of studies have discussed cooling flows as the primary explanation for the observed nuclear 
activity and star formation, through accretion of the cooling filaments into the central galaxy 
(e.g., Cowie et al. 198^ ; Canizares et al. 1987] , Fabian fc Crawford 1990| ). However, cooling flows 



need not be the sole mechanism responsible for triggering and fueling the AGN — in particular, 
the long timescales of a steady cooling flow cannot be reconciled with the comparatively short 
radio source lifetimes without invoking additional mechanisms to trigger the onset of the radio 
source. Furthermore, many AGN with similar properties are found in host galaxies not associated 
with cooling flow clusters. 

Thus we are motivated to discuss the plausibility of an additional mechanism which is 
generally associated with AGN triggering and fueling, namely interactions between galaxies (e.g.. 



Balick &: Heckman 1982 ; Barnes fc Hernquist 1992 , and references therein). Specifically, "minor 



mergers" between the host galaxy and a small gas-rich companion are an important means of 
dumping fuel into the central region ( Hernquist Sz Mihos 199^ ; Walker, Mihos, Sz Hernquist 199^ ; 



De Robertis, Yee, Sz Hayhoe 1998). In such cases the host galaxy retains its original morphology 
more or less intact, with the primary change consisting of a fresh supply of ~ 10^ — 10^ Mq of 
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stars, gas and dust. 

The central elliptical in A2597 displays a number of features characteristic of a recent "minor 
merger": (1) substantial dust which appears dynamically unsettled, located predominantly on 
one side of the nucleus; (2) the extended complex of neutral, ionized and molecular gas, showing 
strong asymmetries on small scales (the present work) and large scales (|Heckman et al. 1989); 



(3) a radio source with an age ^ lO'^ yr; (4) small galaxies nearby, which if dynamically associated 
with the cD would have supplied the gas; (5) the clusters of young stars. The high-pressure cluster 
environment associated with a cooling flow can enhance the star formation to the observed levels 
and confine the radio source. However, the actual triggering mechanism of the star formation 
and AGN is plausibly accounted for by this type of interaction, likely with a galaxy containing 
substantial amounts of molecular gas as discussed in the previous section. When such a galaxy 
interacts with the central elliptical, some of its surviving molecular material would be disrupted 
and become visible as atomic/ionized gas, while the remainder is accreted into the core of the cD. 

In this scenario one expects comparatively little reddening in the diffuse line-emitting 
filaments, due to the relatively short sputtering times of dust in the high-pressure surrounding 
hot gas ( Draine &: Salpeter 1979| ); any remaining dust is concentrated in regions containing 



substantial amounts of dense molecular gas. Indeed, we find that the strongest extinction in A2597 
is concentrated in the compact, dusty regions immediately adjacent to the nucleus, located in the 
same vicinity as the molecular H2 emission detected by Donahue et al. (199^ ). It is therefore 



plausible that the observed distribution and physical properties of the dust and gas in A2597 can 
be accounted for by the capture of an infalling gas-rich galaxy which has retained a significant 
fraction of its original dense clouds of molecular gas and dust. This would also provide a supply of 
fuel for the current phase of activity of the AGN. The age of the radio source is somewhat younger 
than the timescale associated with the interaction, but this likely represents the additional time 
required to funnel gas into the center of the galaxy, which in numerical simulations typically occurs 
well into the merger process; alternatively the currently observed radio source may be only the 
latest in a series of multiple outbursts during the galaxy interaction process. Finally, the location 
of this object in the center of a dense, high-pressure ICM would account for the enhanced star 
formation, as well as the apparent strong confinement of the radio source. 



7. Summary and Final Comments 

In this paper we have presented HST WFPC2 images that spatially resolve some of the most 
intriguing features of the cD elliptical in A2597: 

1. A complex system of emission- line filaments occupies the inner few kpc of the galaxy, 
displaying a number of close spatial associations with the synchrotron-emitting radio lobes. The 
morphologies of several of the filaments are strongly suggestive of direct interactions between the 
radio source and ambient gas. In particular, we find that the detailed physical properties of the 
line-emitting gas around the northern radio lobe can be described in terms of a shell of neutral 
and ionized gas that surrounds the expanding radio lobe, and we show that the energetics of this 
shell can be readily accounted for by energy input from the expanding radio lobe. We also resolve 
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a knot of high surface-brightness Une emission approximately co-incident with a radio hotspot 
to the south-west of the nucleus, and we show that its properties are consistent with those of a 
dense, inhomogeneous clump of gas (nj-ms ^ 14 cm~^, M^ ^ 6 x IO^Mq) sufficiently massive to 
deflect the radio jet over the lifetime of the radio source and account for the severely disrupted 
morphology of the southern radio lobe. Finally, the lack of detectable line emission in a number of 
locations around the radio lobes suggests an intrinsically asymmetric distribution of ambient gas 
about the galaxy. 

2. The excess blue continuum emission is clearly resolved by our HST images into a number 
of components, specifically an intense band of blue emission associated with the base of the 
northern radio lobe, as well as blue clusters south-west from the southern radio lobe, and a 
number of more diffuse regions distributed about the radio lobes. The morphological properties of 



the blue continuum (together with its low observed polarization, McNamara et al. 1998| ) rule out 
several explanations for its origin, including optical synchrotron or inverse Compton emission, or 
scattered nuclear light, leaving recent star- forming regions as the most likely remaining possibility. 
Comparing our photometry with stellar population synthesis models shows that these regions 
can be modeled either as single-burst populations with ages of several times 10^ — 10*^ yr, or as 
continuous (constant) star formation with ages about an order of magnitude higher, depending 
upon details of the assumed metallicity and IMF parameters. 

3. The strong dust obscuration in the central region is resolved into a number of large 
clumps, each several hundred pc in size. Most of these are distributed in a region extending away 
from the nucleus toward the south-east, while another strong clump of obscuration is evident 
immediately to the north of the nucleus, and some less heavily absorbed regions are present just 
west of the bright radio / emission-line hotspot in the south. A striking feature of the dust is its 
strong degree of clumping and highly asymmetric spatial distribution about the galaxy center. 

The cD elliptical in the cooling flow cluster A2597 is a complex system: it contains a 
powerful, compact radio source, surrounded by an extensive network of emission-line fllaments 
with significant amounts of dust and prominent regions of recent or on-going star formation. We 
suggest that the high-pressure environment associated with the cooling flow is crucial to enhancing 
the star formation rate to the high levels observed. However, while a steady cooling flow can also 
provide a source of fuel for the star formation and the AGN, its long timescales are difficult to 
reconcile with the short radio source lifetime, the substantial amounts of dust and molecular gas 
in the central few kpc, and the highly asymmetric distribution of the gas, dust and young stars 
without invoking an additional, more recent mechanism to account for these features. 

We propose that a recent "minor merger" or interaction in which the cD accreted gas and dust 
from a smaller gas-rich galaxy provides a plausible mechanism, specifically in terms of triggering 
the formation of the radio source and introducing large amounts of molecular gas and dust into 
the cD elliptical. The observed asymmetry of the gas and dust would indicate that the system 
is still in the process of virialization, i.e., the event took place within the last ~ 10^ — 10^ yr. 
This dynamical timescale is also compatible with the timescales derived from the stellar models 
involving continuous star formation with a constant SFR, i.e., a scenario in which the interaction 
is responsible for inducing on- going star formation. However, if the star formation is described 
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in terms of a single burst then the models imply ages of the order ^ lO'^ — 10^ yr, in which case 
the star formation would more likely have been triggered by the passage of the radio source. 
The age of the radio source is substantially less than the galaxy interaction timescale, and this 
either represents the additional time required to funnel gas into the center of the cD, or otherwise 
indicates that the present phase of radio activity is the latest in a series of multiple outbursts. 

The properties of A2597 are directly tied to the general question of the relationship between 
galaxy interactions, nuclear activity and star formation. If active galaxies in cooling flows such 
as A2597 are triggered in the same way as AGN unassociated with cooling flows, then the 
primary reason for the difference in their properties, specifically the compact radio morphology 
and enhanced star formation, results from being located at the center of a dense, high-pressure 
ICM. Better understanding the effect of environment upon star formation would be especially 
desirable in the context of studying objects in a range of environments at low redshifts, as well as 
investigating the formation of stars and galaxies at early epochs where the environment is likely 
to play an important, perhaps dominant, role. 
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Fig. 1.— Broad-band F702W HST WFPC2 Planetary Camera image of A2597. The spatial scale 
is ~ (y.'0455/pixel, corresponding to ~ 65 pc / pixel for our adopted cosmological parameters. The 
orientation of the image is such that north is to the top and east is to the left, and the original 
pointing orientation of the image has been preserved (in order to avoid smoothing the image by 
rotating the pixels). 
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A2597 - Greyscale: HST/WFPC2 [0II]A3727; Contours: VLA 3cm 
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Fig. 2.— Narrow-band F410M HST WFPC2 WF3 image of A2597 (greyscale). The spatial scale 
is ~ O'.'0996/pixel, corresponding to ~ 140 pc / pixel for our adopted cosmological parameters. At 
the redshift of the object, this filter fully samples the [O II] A3727 emission line, and this image 
can be considered representative of the emission-line structure, although we discuss in the text 
the possibility that continuum emission also contributes to some of the filament structure. We 
superpose contours from a 3.5 cm VLA radio continuum image obtained by Sarazin et al. (1995| ). 
Note the extremely good morphological agreement between the northern radio lobe and emission- 
line arc, and also between the southern radio hotspot and the bright, compact line-emitting knot. 
The southern emission-line filament also lies near the eastern edge of the southern radio lobe. 
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Fig. 3. — Broad-band ?7iF450W — ''7^F702W color distribution of the galaxy, obtained from the PCI 
F702W and F450W images, each of which was convolved to the resolution of the other before 
converting to magnitudes and performing the subtraction. The magnitude units have been retained 
in the HST VEGAMAG system. We superpose contours from the 3.5 cm VLA radio continuum 
image ( parazin et al. 1995| ). 
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Fig. 4. — Stellar evolution models (Chariot & Bruzual 1999), K-corrected to the redshift of A2597 
and calibrated for our HST WFPC2 filter bandpasses. In the top panel we plot the color evolution 
of models with constant SFR, for Salpeter and Scalo IMFs, as well as instantaneous burst models 
with a Salpeter IMF and metallicities of Zq, 0.2Zq and 0.02Zq. For purposes of clarity, the 
instantaneous burst models with a Scalo IMF are omitted from the plot since their color evolution 
is very similar to the Salpeter models. The constant SFR models, on the other hand, show little 
difference between Zq and 0.2^0 models, therefore we plot only 0.2Zq and 0.02Zq for these. 
The lower panel displays the evolution in color-magnitude space for instantaneous burst models of 
10^ Mq with a Salpeter IMF and metallicities of Zq, 0.2Zq and 0.02Zq (represented by the same 
line styles as in the top panel). Again, Scalo models are omitted, the color differences between the 
Scalo and Salpeter models being much less than the formal uncertainties on our photometry. We 
overplot the measured colors and apparent magnitudes of the four blue compact objects south-west 
from the southern radio lobe. The ages of the model bursts are indicated on the plot, and the 
inferred ages of the clusters are discussed in the text. 
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Fig. 5. — The final set of isophotal parameters fitted to tlie galaxy continuum distribution, together 
with the mF450W — "T'F702W color distributions. Error bars are absent from the position angle and 
ellipticity plots since these parameters were kept constant in producing the final fits to the F702W 
and F450W images. 
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Fig. 6. — The greyscale represents the "residual" filamentary structure in the F702W band, 
obtained by subtracting the elliptical galaxy model from the original image. The contours are 



the 3.5 cm VLA radio continuum image ( Sarazin et al. 1995 ) 
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